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Abstract 

We present results from Molecular Dynamics simulations of the thermal glass 
transition in a dense polymer melt. In previous work we compared the simulation 
data with the idealized version of mode coupling theory (MCT) and found that the 
theory provides a good description of the dynamics above the critical temperature. 
In order to investigate the influence of different thermodynamic paths on the struc- 
tural relaxation (a-process), we performed simulations for three different pressures 
and are thus able to give a sketch of the critical line of MCT in the pressure- 
temperature-plane [(p, T)-plane] , where, according to the idealised version of MCT, 
an ergodic-nonergodic transition should occur. Furthermore, by cooling our system 
along two different paths (an isobar and an isochor), with the same impact point 
on the critical line, we demonstrate that neither the critical temperature nor the 
exponent 7 depend on the chosen path. 

PACS: 61.20.Ja,64.70.Pf,61.25.Hq,83.10.Nn 
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1 Introduction 

The understanding of the glass transition has been a longstanding problem of condensed 
matter physics and materials science [III [2, [5[ f|]. The term glass transition is used to 
describe a phenomenon, where a solidification of the liquid without simultaneous crys- 
tallisation occurs. In a narrow temperature region the viscosity of the material increases 
by some 14 orders of magnitude, and, although the typical time scales of the system 



change tremendously, no significant change in the structure is observed and the material 
remains amorphous. 

Experimental research in this area of scientific interest has been conducted for more 
than a hundred years |], || . Since then a large number of phenomenological theories were 
proposed, such as the free volume theory |6|, [7], |8| or the Gibbs-Di-Marzio theory || [II], 



PI , |P2]1 , which tried to explain the characteristic features of the glass transition. These 
theories developed models, mainly based on physical intuition, from which thermodynamic 
properties and the temperature dependence of the viscosity, i.e., the Vogel-Fulcher law 
IQ, could be derived. However, an explicit relationship between the model parameters and 
the microscopic properties of the glass former remained hard to establish, and a detailed 
description of the shape of dynamic correlation functions was not attempted. 

In recent years the mode coupling theory (MCT) was successful in describing a broad 



range of features observed in experiments [13L O, OH [Tq, 17] and simulations [18, 1£. 20, 



2T , |22|, 23, |23|, ^3], (also see [27J for an overview) in a temperature region close to a 
critical temperature T c , which in general lies above the empirically defined glass transition 
temperature T g . The critical temperature is associated with a change in dynamics from 
a liquid-like to a solid-like behaviour without any concomitant modification of the glass 
former's structure. In contrast to most earlier theories this theory starts from well known 
microscopic dynamics and uses techniques already applied in the field of critical dynamics 
to derive a set of dynamic equations for the density correlation functions of the system 
p8| , |29| , |5C\ pTf (also see |32|, |3j| for review articles). It made novel predictions on the 
temporal behaviour of the density correlations, which could be tested experimentally and 
generated new research in this 

The ideal version of the theory states that at the critical temperature the glass former 
should show a transition from ergodicity to nonergodicity P5| |. However, in real systems 
this transition is in general absent (an exception being (probably) some colloidal glasses 
|]T4l|), because additional relaxation processes, which are not included in ideal MCT, assure 
that the system is able to relax even below the ideal transition temperature. In contrast, 



the extended version of MCT [29, 30, 31] does account for these additional processes in 



terms of activated hopping processes and therefore predicts that the transition is indeed 
avoided. Since sufficiently above the critical temperature the contribution of jump pro- 
cesses to the dynamics of the glass former is negligible, the ideal version of MCT often 
describes experimental data in this temperature regime very well. Only if one is suffi- 
ciently close to or even below the critical temperature, one has to apply extended MCT 



13, 19]. However, such temperatures are seldomly accessible in computer simulations of 



polymer systems, at least if equilibrated melts are required. Thus, in the following section 
we will concentrate on the predictions of ideal MCT. 

Although MCT has been applied to experimental data numerous times, many aspects 
of the theory still remain to be thoroughly investigated. Little is known on the influence of 
external thermodynamic parameters on the transition temperature (although see [34. 35. 
3"Bl |37J) or the a-process. Since in general higher pressure causes higher densities, which in 



turn means that the movement of an individual particle is more hindered, an increase in 
pressure results in an increase of the transition temperature. This effect is well known for 
the glass transition temperature ||. Hence, it is possible to cause a glass transition by sole 
increase of the pressure, which indeed has been observed in experiment p7| . In general 
the thermodynamic parameter space (for instance, the (p, T)-plane) should decompose 



into two fluid phase and an ideal glass phase, both separated by a critical line 

at which the transition to nonergodicity should occur. Since MCT applies to systems in 
thermodynamic equilibrium, the position of the critical line should be independent on 
the thermodynamic path chosen upon cooling. Furthermore, close to the critical line, the 
exponents of the theory, which determine most of the quantitative behaviour of the glass 
former, depend solely on the impact point of the thermodynamic path on the critical line, 
granted one does not choose a too exotic path, e.g. one that runs almost parallel to the 
critical line. Therefore two different thermodynamic paths, which have the same impact 
point on the critical line, should not only yield the same critical temperature, but the 
quantitative behaviour of the systems, described in terms of MCT, should be the same 
along both of them. To our knowledge, so far this prediction has never been verified. 

Thus, we decided to conduct a study of the influence of pressure on the parameters 
of ideal MCT, which also should enable us to give a sketch of the critical line in the 
(p, T)-plane. To this end, we chose a model for a polymer melt, which has already been 
used in an earlier study of the glass transition [p6fl . Clearly, because of the connectivity of 
the monomers along the chain, our model is by no means a simple liquid, and therefore 
certainly not the kind of system ideal MCT was originally developed for. On the other 
hand, polymers, according to Angell's classification scheme J38|, mostly belong to the 
fragile glass formers, to which MCT has been applied successfully, and are extremely 
good glass formers, i.e. "supercooled" melts in thermal equilibrium can be prepared very 
well. In previous work |26| it was demonstrated that it is possible to equilibrate our model 
well in the regime of the supercooled melt, a prerequisite to apply MCT, and that the a- 
relaxation behaviour is compatible with MCT. Since we were able to simulate the system 
in the isochoric (NVT) as well as in the isobaric (NpT) ensemble, we try to test the 
prediction of thermodynamic path independence in the present paper. To this end, we 
cooled our system along an isochoric path, which shared the impact point on the critical 
line with one of the isobaric paths. As in our earlier work ||26fl , we concentrate on the 
a-relaxation behaviour in this study. An analysis of the /3-relaxation can be found in 



Ref. gg]. 

The remainder of the paper is organised as follows. In section II we briefly describe our 
model and the simulation technique. By performing simulations along different isobars we 
were able to investigate the influence of pressure on the critical temperature. The results 
of these simulations will be discussed in section III. In section IV we will perform a test 
of the thermodynamic path independence, and in section V conclusions will be drawn. 



2 Model and Simulation Technique 



For modelling the inter- and intramolecular forces we used a bead-spring model derived 



from the one suggested by Kremer and Grest |^0| and also used in several recent simu- 
lations El However, we included here also the attractive part of the Lennard- Jones 



potential, since previous work on a lattice model for a glassy polymer melt p3| , [44|| had 
shown that without such an attraction the model would produce a negative thermal ex- 
pansion coefficient. The model of Kremer et al. EO, fill E3] is close to an athermal model 



of polymer melts and hence does not exhibit a glass transition driven by temperature at 
all. 



As in our past simulations each chain consisted of 10 beads with mass m set to unity. 
Although these chains are rather short, they already show the static behaviour charac- 
teristic of long polymers in the melt (e.g. Gaussian statistics for the end-to-end distance 
distribution, a Debeye scattering law for the single chain structure factor, etc.). Note that 
each bond in this model would correspond to n « 3 — 6 covalent bonds along the backbone 
of a real chain, if one were to map this coarse-grained model onto a real polymer. Between 
all monomers there acted a truncated Lennard- Jones potential: 
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where C was a constant which guaranteed that the potential was continuous everywhere. 
Since it was not our aim to simulate a specific polymer, we used Lennard- Jones units, 
where e and a are set to unity. Note that this means that all quantities are dimensionless. 
In addition to the Lennard- Jones potential a FENE backbone potential was applied along 
the chain: 
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The parameters of the potential were taken as k = 30 and Rq = 1.5, guaranteeing a 
certain stiffness of the bonds while avoiding high frequency modes and chain crossing. 
Furthermore, with these parameters we set the favoured bond length to a value slightly 
smaller than the length favoured by the Lennard- Jones potential. Thus we introduced two 
different incompatible length scales in our system, which prevents the emergence of long 
range order (i.e. crystal formation) at lower temperatures. 

Unlike previous lattice models for the thermally driven glass transition of polymers 



43] , [44]], the present model has a qualitatively reasonable equation of state with a pos- 
itive thermal expansion coefficient, and can easily be studied under constant density or 
constant pressure. It allows to study motion and structure from local scales (motions in 
the neighbour cage) upto large scales. 

In order to keep the temperature fixed, all simulations were performed using a Nose- 



Hoover thermostat [^, IB, 4"7fl . In this technique the model system is coupled to a heat 
bath, which represents an additional degree of freedom. To set the system to a desired 
pressure, the size of the simulation box was adjusted to yield the correct density at each 
temperature. The resulting configurations were used as start configurations for runs in the 
canonical ensemble, where the size of the simulation box was kept fixed. Only during these 
canonical runs dynamic correlation functions for the further analysis were calculated. A 
more thorough discussion of the simulation technique applied can be found elsewhere 



26, 58 1 . Here we only emphasize that we have carefully checked that the Nose- Hoover 
thermostat does not lead to any artefacts in the dynamics of the single chain correlators 
and local properties that were studied here [f4"8"|. Note also that our chain length iV = 10 



was short enough that our results are not affected at all by chain entanglement effects. 

Altogether we performed simulations at more than 40 different points in the thermo- 
dynamic phase space. At each point 10 independent configurations were simulated, each 
consisting of 120 polymer chains of 10 monomers. With this we were able to perform 



simulations along three isobars and one isochor. Table 1 shows which temperatures were 
simulated in which ensemble. 

In order to equilibrate an individual system at lower temperatures one had to simulate 
for very long times (> 10 6 MD-steps). Generally, the equilibration of the lowest temper- 
ature in a given ensemble lasted as long as the sum of all equilibration times at higher 
temperatures of the same ensemble. Altogether the simulations consumed an equivalent 
of approximately 10 CPU-years on a PentiumPro processor run at 180 Mhz. 



3 Dynamical properties at different pressures 

As already discussed in the opening paragraph the glass transition manifests itself by a 
steep increase of the relaxation times by several orders of magnitude. In order to extract 
these time scales from the simulation data we computed a number of dynamical quantities, 
like the incoherent intermediate dynamic structure factor: 

W = (^E^ (r,(t) " r4(0)) ). ( 3 ) 

where N stands for the total number of monomers in the melt. This function measures the 
self correlation of the particle positions at different times, and, by varying the wave- vector 
q, at different length scales. 

Recently, orientational degrees of freedom and their relaxational behaviour have come 
into focus of theoretical research on the glass transition [^, 5(|. Results of molecular 



dynamic simulations for a fluid consisting of diatomic molecules |22], ^3], 23J] illustrated 
that there can be significant differences between orientational and translational relaxation. 
Such differences are also observed in experiments (see Ref. for instance). Clearly, it 
should be interesting to check, whether we could find any differences between orientational 
and translational relaxation in our model. Hence, we also calculated the orientational 
correlation of the end-to-end vector: 

sWlJ ,, It 6 '"' „ )), ,= 1,2,..., (4) 



e(t) mi e(0) 

where L n stands for the nth Legendre-polynomial, e(t) is the end-to-end vector of a 
polymer at time t, and || e || is the length of the end-to-end vector at time t. The same 
formula can be applied to measure the dynamical correlation of a bond vector h(t): 

Equations @ and (f|) characterize the reorientation dynamics of the largest and of the 
smallest vectors along the backbone of a chain. In the analysis, we only calculated the first 
and second polynomial, since these quantities can be measured by dielectric relaxation 
and light scattering, respectively. 

With these three dynamical correlation functions we define the following correlation 
times: 

^(r 9 ) = 0.3 E n {r En )= 0.3 B n (r Bn ) = 0.3. (6) 



We have computed a number of other related quantities as well, such as the Rouse-modes 
of the system or the mean-square-displacements, which are discussed in other publications 

© H 0- 



3.1 Behaviour of dynamical correlators 

As shown in Fig. [I] and Fig. [2], the correlators decay in one step at high temperatures, 
while at lower temperatures a two step process starts to emerge, which becomes the 
more pronounced, the lower the temperature. The emergence of a plateau in the decay is 



related to the cage effect , where an individual monomer is trapped by its surrounding 
particles. The average time a monomer needs to escape from the cage of its neighbours 
increases with decreasing temperature, which explains the extension of the plateau. The 
presence of a two step relaxation, the so-called /3-relaxation (onto and off of the plateau) 
and a-relaxation (off of the plateau and long-time structural relaxation), is a common 
feature of glass formers, and is also predicted by MCT. 

As can be seen from the plots, the qualitative behaviour is not affected by the applied 
pressure, although at higher pressure the two step process starts to show up at higher 
temperatures already. Furthermore, while the height of the plateau depends on the specific 
correlator, it hardly varies with pressure. The orientational correlators [first and second 
Legendre-polynomial, see Eqs. (|J) and @] exhibit a rather high plateau value which is 
often close to unity (the plateau of the first Legendre-polynomial is always larger than that 
of the second). Therefore the two step process is only visible on magnification. Clearly, 
the contribution of the a-process to the overall relaxation of a correlator depends on the 
quantity considered. 

Another characteristic of glass forming liquids is that close to the critical tempera- 
ture the time temperature superposition principle should hold for the a-relaxation. One 
therefore has to rescale a dynamical correlator by a suitably defined a-relaxation time 
and to check whether the curves fall on a master curve in the a-regime. As we reported 
in our earlier work for p = 1 (and constant volume) , this is indeed the case. Here, we 
additionally observe that our data also obey a time-temperature-pressure superposition 
principle, i.e., in the a-regime data taken from different isobars collapse on a single master 
curve. This is illustrated for a number of different dynamic correlators in Fig. RL where 
fourteen different curves are included in one plot. Similar behaviour has been observed 



in experiments on orthoterphenyl P5 |, but, to the best of our knowledge, this is the first 
report from computer simulations for such a behavior. 



3.2 Behaviour of relaxation times 

Figure |I| and ^ show that on lowering the temperature an increase of the relaxation 
times by several orders of magnitude takes place, as expected for a glass forming liquid. 
By means of asymptotic expansions the idealised MCT derives a number of predictions 
concerning the behaviour of the liquid close to the critical temperature. One finds that 
sufficiently close to T c the increase of the a-relaxation times can be described by the 
following formula: 

t = t°(T-T c )-\ (7) 



where r° is an amplitude which depends on the specific relaxation time considered, and 
7 is a parameter of the theory which should be the same for all correlation times, if the 
corresponding correlator couples to density fluctuations. Furthermore our analysis of the 
/^-regime suggested that 7 should take the value of 7 = 2.09 for the isobar p = 1.0 |39| . 

At all pressures investigated, it is indeed possible to locate a temperature interval, 
where the increase of the relaxation times can be described by Eq. ([7|). When applying 
Eq. (|7|), r°, T c , and 7 were treated as adjustable parameters. Although it is not incom- 
patible with MCT that rotational and translational degrees of freedom freeze at different 
state points in the temperature-density plane, as was demonstrated in recent publica- 



tions |49|, |50[ , our analysis suggests that it is possible to find a critical temperature for 
all isobaric paths, which is independent of the specific correlator and solely a function of 
the pressure considered. The critical temperatures and densities we obtained are listed in 
Table |. Note that the error for the critical temperature for p = 1.0 is smaller than for the 
other pressures because a larger number of temperatures was simulated for this isobar. 
The pressure dependence of the critical temperature is depicted in Fig. £2[ As expected, 
the critical temperature increases with increasing pressure as also calculated for Lennard- 
Jones models in [|34| . As one can also see from Table [| within the error bars the quantity 
p c T~ 1//4 is a constant at the mode coupling critical point as also found experimentally in 



e.g. [37] and in the simulation of soft sphere models [51, 64 53] and for Lennard- Jones 



mixtures 54 . The value we found is within the error bars identical to the Lennard- Jones 



value in [B3 



Figure y shows a double logarithmic plot of a-relaxation times against T — T c , using 
the critical temperatures of Eq. (^). For all pressures there is a temperature interval, 
where the data points lie on a straight line in accord with Eq. (|7|). Deviations from 
the power-law behavior are visible both at small and large distances from the critical 
point. The deviations for large T — T c are expected because Eq. ([fj) is an asymptotic 
expansion which is only valid, if the reduced distance to T c , i.e., (T — T c )/T c , is small. 
Upper bounds for the validity of Eq. (|7|) are approximately 0.7 (p = 0.5), 1.2 {p = 1), 
and 0.6 (p = 2), which is comparable to experiments |13|, [15! anc ^ other simulations 
fl2D| , PL |23, However, these upper bounds strongly depend on the quantity under 

consideration. Whereas deviations are very pronounced for the smallest length scale (q = 
9.5), Eq. (|7|) provides a good description at all, except perhaps at the lowest temperatures 
for the end-to-end distance. 

On the other hand, the deviations from the idealised power law at low temperatures 
could be attributed to the ergodicity restoring jump processes mentioned above. Close to 
the critical temperature these processes start to contribute significantly to the relaxation 
dynamics of the system, and therefore the actual relaxation times can be smaller than 
the predictions of idealised MCT. This behavior has been discussed in experiments, e.g. 
[p~3] , p~5|[ , and simulations, e.g. |Lj], |22], p5fl . Therefore, in practical application of Eq. ([?[) 



one faces the problem that its range of validity is limited from below and above, and that 
it additionally depends on the quantity under consideration. 

Furthermore, Eq. (|7]) implies that in the temperature regime, where the idealised 
MCT is applicable, the ratio of two different a-relaxation times should be independent of 
temperature. As demonstrated in Fig. |6], this is not the case, even in the regime where the 
/3-analysis could be done, i.e., for T — T c < 0.07 | 39|| . The ratio between different relaxation 
times can change by almost a factor of two, and the effect is stronger for q = 9.5 than for 



q = 6.9 (first minimum and maximum of the static structure factor, respectively). Note 
that we obtain the same result when applying a different definition of the a-relaxation 
time, which includes the nonergodicity parameter f q c , i.e., 4> q (r q ) = e _1 /g C - However, 
it is not clear whether this finding is a strong contradiction to MCT, because we have 
eliminated the dominant temperature dependence, given by Eq. (|7|), when dividing two 
relaxation times. Since we are close, but not very close to T c , and Eq. (0) is, strictly 
speaking, only asymptotically valid, one could expect a smooth temperature dependence 
of prefactors. Such a conclusion can also be drawn from Ref. | T6fl , in which the MCT- 



equations for a model of a colloidal suspension are solved numerically and compared with 
the asymptotic results. There, it is found that the ratio of two relaxation times becomes 
constant only very close to the critical point, although Eq. (0) is already observed for 
larger distance to the critical volume fraction (see Fig. 7 of Ref. |16[). Interestingly, Fig. | 
shows that the ratio is not a monotonic function and exhibits a maximum approximately 
at the beginning of the temperature interval, in which we can apply ideal MCT to describe 
the a-relaxation time. It seems as if at this temperature a change in the dynamics of the 
system occurs. 

This problem is also reflected in Fig. [7] which shows the results for 7 when fitting 
Eq. (0) to the a-relaxation time of <p s q (t) at p = 1. The critical temperature was kept fixed 
(T c = 0.45) in the fits, and the maximum possible number of temperatures was taken into 
account to determine 7. It is interesting to note that the 7- values, determined from Fig. [| 
for the different pressures, agree with one another within the error margins so that the 
following discussion is not specific for p = 1. Figure [7] shows that the fit procedure yields 
a decrease of 7 with decreasing q, but the 7-values are distributed around the result of 
the /3-relaxation analysis, 7 = 2.09 p9| . Alternatively one can keep the exponent 7 = 2.09 
constant and adjust the critical temperature, T c m, p9|| . Then the critical temperatures for 
q > 3 coincide within the error bars with the value obtained from the /3-analysis. However, 
the diffusion coefficient of a chain yields a T c , which is significantly lower. Physically, both 
types of analysis suggest that going from T c to higher temperatures the melt has a stronger 
tendency to liquify on short than on the large length scales. Such a behavior is not unique 
to our polymer model, but It was also found in other simulations ET 



4 Test of the thermodynamic path independence 

In order to verify the prediction of the thermodynamic path independence, we estimated 
the density of the melt at the critical temperature (T c = 0.45) of the isobar p — 1. The 
density is p = 1.042. Then we performed a number of simulations in the NVT-ensemble at 
the appropriate isochor (schematically, this is illustrated in Fig. |j), and again calculated 
various dynamic correlation functions at the simulated temperatures (see Table [[]). 

The qualitative behaviour of the dynamical correlation functions along this isochor 
does not differ from the behaviour observed at the various isobars, discussed in the pre- 
ceding section. As can be seen in Fig. ^| and Fig. pi], which show exemplarily the dynamic 
correlation of the end-to-end vector orientation, we find again that at lower temperatures 
a two step relaxation occurs (which cannot be seen on the scale of figure due to the large 
plateau value), that the relaxation times show a steep increase, and that at least for the 
lower temperatures the time temperature superposition principle holds. This could have 



been expected, since also earlier simulations of the model in the NVT ensemble had shown 



such a behaviour |ftq| . It is interesting to note, however, that, compared to the appropri- 
ate isobar, the two step relaxation process can now be observed at higher temperatures 
already and that for the studied temperatures the a-relaxation time is almost one order 
of magnitude larger. 

Once again it is possible to find a temperature region, where the behaviour of the 
a-relaxation times, extracted from the different dynamic correlation functions, can be 



described by Eq. (0). This is illustrated in Fig. [□], where we show the temperature de- 
pendence of various correlation times, as measured in the NVT-ensemble, plotted in such 
a way that the applicability of the MCT-prediction is clearly demonstrated. Qualitatively, 
we find the same features, as discussed before for the isobars. There are deviations from 
linearity at large temperatures, the deviations are more pronounced for the shortest length 
scales, but still the fits yield very similar values for the critical temperature, which can 
be combined to: 

T c (p= 1.042) = 0.445 ±0.010 . (8) 

Within our error bars the value coincides with the critical temperature obtained for the 
isobar p = 1.0. For the NVT-simulation the error bar is larger, since we were not able 
to equilibrate the melt as closely to the critical temperature as it was possible in the 
NpT-simulation. Note that the lowest temperature in Fig. |5|b is T = 0.46, whereas it is 



T = 0.5 in Fig. [IT]. This difference is caused by the larger relaxation time in the NVT- 
ensemble (due to higher density/pressure) at a specific temperature in comparison to the 
NpT-ensemble. Therefore the estimate becomes less accurate, but we can still conclude 
that the critical temperature of ideal MCT is indeed independent of the thermodynamic 
path chosen. 

Finally, we want to verify whether the exponent 7 is independent of the thermodynamic 
path as well. As already discussed in section III, the exponent 7 shows a pronounced 
dependence on the dynamic correlation function considered, if one works with the same 
critical temperature for all quantities and extends the fit interval as much as possible. The 
same dependence is also found for the isochor, but the results coincide within the error 
bars with those of the isobars, as Fig. |TT] illustrates. Therefore, 7 is in fact independent 
of the chosen path, which exemplifies the thermodynamic character of the critical point 
in mode-coupling theory. 



5 Conclusions 

In this paper we have presented results of a large scale molecular dynamics simulation 
for a supercooled polymer melt. Our model is a coarse-grained bead-spring model with 
nonlinear springs connecting monomers along a chain and Lennard- Jones interactions 
between all monomers. By including competing length scales in the model we prevented 
the melt from crystallising at lower temperatures. 

The present study concentrated on the influence of pressure on the a-relaxation be- 
haviour of the melt. Upon cooling we see a steep increase of the a-relaxation time, 
and all dynamic correlation functions show a two-step relaxation. By comparing data 



from different isobars, we found that our system does not only exhibit time-temperature- 
superposition above T c , but time-temperature-pressure-superposition as well. 

For all pressures investigated it has been possible to locate a temperature interval, 
where the increase of the a-relaxation times could be described by idealised MCT. There- 
fore we have been able to investigate the dependence of the critical temperature of MCT 
on pressure and to give a sketch of the critical line in the (p, T)-plane. However, whereas 
the critical temperatures, determined from different quantities probing both small and 
large length scales of the melt, coincide within the error bars, the approach towards T c , 
i.e., the exponent 7, is very sensitive to the precise choice of T c in the fit, and depends 
on the quantity considered. When fixing T c , we find that the a-relaxation times of (jf q (t) 
for g-values distributed around the maximum of the structure factor are compatible with 
result of the /^-analysis [p9| . Deviations occur for much smaller and larger q- values. The 
deviations at large q can be explained by the sensitivity of 7 on T c , since fixing 7 at the 
value of the /3-analysis instead of T c yields estimates for T c that are compatible with the 
result of the /^-analysis p9 |. However, such an alternative fit procedure does not remove 



the discrepancies found on the largest length scale. On these length scales, 7 is smaller 
than expected from the /3-analysis. Similar deviations are also observed on smaller length 
scales, if the critical point is approached very closely. They can be rationalized, within 
the theoretical framework of MCT, by ergodicity restoring processes, which compete with 
and finally dominate over the cage effect treated by the idealised theory, if T < T c . To 
what extent the predictions of the idealised theory are observable, therefore depends not 
only on the quantity under consideration, which is also pointed out in recent theoretical 
work [|T^, [57|], but also on the distance to the critical point. If one is too close, ergodicity 
restoring processes interfere, and if the temperature is too large, one leaves the asymptotic 
regime, where the formulas of the idealised MCT are expected to hold. 

By performing simulations along an isochor which had the same impact point on the 
critical line as one of the isobars, we have been able to verify, that within the error margin 
the critical temperature of MCT is indeed independent of the thermodynamic path one 
chooses upon cooling. Furthermore, we have shown, that the exponent 7 does not depend 
on the thermodynamic path either, within the caveats explained in the last paragraph. 

In summary, one can therefore say that the idealised theory is a good starting point for 
a quantitative description of the dynamics above T c , and seems to capture the essential 
physics, not only for simple liquids, but also for polymer model. Reasons, why this could 



be the case, are further discussed in Ref. 39 
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Table 1: The table shows, at which temperatures and densities or pressures simulations 
were performed. 



ensemble 


simulated temperatures 


isochor (p = 1.042) 


0.5, 0.52, 0.55, 0.58, 0.6, 0.65, 0.7, 0.8, 0.9, 1.0, 2.0 


isobar (p = 0.5) 


0.45, 0.48, 0.5, 0.52, 0.55, 0.6, 0.7, 1.0 


isobar {p = 1.0) 


0.46, 0.47, 0.48, 0.49, 0.5, 0.52, 0.55, 0.6, 0.65, 0.7, 1.0, 2.0, 4.0 


isobar (p = 2.0) 


0.52, 0.55, 0.57, 0.6, 0.7, 0.8, 0.9, 1.0, 2.0 



Table 2: Critical temperatures and densities and soft sphere scaling variable at the critical 
point 



p 


T 

- 1 c 


Pc 




0.5 


0.425 ±0.010 


1.035 ±0.01 


1.28 ±0.02 


1.0 


0.450 ±0.005 


1.042 ±0.01 


1.27 ±0.02 


2.0 


0.490 ±0.010 


1.054 ±0.01 


1.26 ±0.02 
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Figure 1: Intermediate dynamic structure factors at the first maximum of the static struc- 
ture factor (q = 6.9) |3{| measured along the isobars p = 0.5 (a), p = 1.0 (b) and p = 2.0 
(c). The broken line shows the value, which we used to define the a-relaxation time scale. 
From right to left, the temperatures decrease, as specified in Table |l|. 




Figure 2: Dynamic correlation functions of the orientation of the bond vectors [second 
Legendre polynomial, see Eq. (|5])]. Different pressures are shown: p = 0.5 (a), p = 1.0 (b) 
and p = 2.0 (c). From right to left, the temperatures decrease, as specified in Table [TJ. 




Figure 3: Figure. 3a is a compilation of the results from Figs, la-c, but with times scaled 
by the a-relaxation time scale for selected temperatures close to the critical temperature 
at the respective pressure: T = 0.45,0.48,0.5,0.52,0.55 for p = 0.5 (dashed lines), T = 
0.46,0.47,0.48,0.5,0.52 for p = 1 (solid lines) and T = 0.52,0.55,0.57,0.6 for p = 2 
(dotted lines). In the a-regime the curves for different temperatures and pressures all 
collapse on a single master curve, demonstrating time-temperature-pressure superposition. 
Figures 3b and 3c show the same behavior for the orientational correlation functions B2(t) 
and E 2 (t) (second Legendre polynomial) of the bonds and the end-to-end vector. Note 
that the plateau for E 2 (t) is so close to 1 that the first step cannot be seen on the scale 
of the figure. 
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Figure 4: Mode Coupling critical temperatures at different pressures. The critical tem- 
peratures represent averages which are derived by fitting Eq. (0) to all relaxation times 
shown in Fig. [| . The broken line is an illustration of the critical line of MCT (guide to 
the eye only), while the arrow symbolises a thermodynamic path at constant density. 
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Figure 5: Temperature behaviour of different relaxation times, measured along the isobars 
p = 0.5 (a), p = 1.0 (b) and p = 2.0 (c). In the plots r q and te 2 are the a-relaxation 
times of the incoherent dynamic structure factor at different wave numbers and the dy- 
namic orientational correlation of the end-to-end vector (second Legendre polynomial), 
respectively. The values of T c are taken from Eq. (0). The solid lines are power-law fits, 
including the largest possible number of temperatures. For p = 1, the fit for q = 6.9 uses 
7 = 2.09, i.e., the 7- value resulting from an analysis of the /3-relaxation . 
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Figure 6: Ratio of different a-relaxation times as measured along the isobar p = 1.0. As 
can be seen, even close to the critical temperature (T c = 0.45) the ratio changes by almost 
a factor of two. 
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Figure 7: Variation of 7 (p — 1.0) with the magnitude q of the wave- vector when fitting 
the a-relaxation time of (jf q (t) by Eq. (0) while keeping the critical temperature fixed 
(T c = 0.45). 
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Figure 8: Dynamic correlation of the end-to-end vector orientation (second Legendre poly- 
nomial), as measured at constant density along the thermodynamic path p = 1.042. Tem- 
perature decreases from right to left, as specified in Table |l|. 




Figure 9: a-scaling plot of the end-to-end vector correlation function for temperatures 
ranging from T = 0.5 to T = 0.7 (see Table |l] for details). The simulation data are for the 
same isochor as in Fig. [S], and the relaxation time was determined by Eq. (|B|). 
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Figure 10: Temperature behaviour of different a-relaxation times, as measured in the 
NVT-ensemble. The solid line represents a fit with Eq. (0). 
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Figure 11: Values of 7, determined from the temperature dependence of various correlation 
times, for two different thermodynamic paths, which yield the same critical temperature. 
On the abscissa, the relaxation times are quoted, from which 7 was determined. For both 
the isobaric and the isochor path the error margins are about 10%, which is rather large, 
since 7 is very sensitive to a variation of the critical temperature. Within the error bars, 
7 does not depend on the thermodynamic path. 



